Temperature Dependence of the Diffusive Conductivity for Bilayer Graphene 



O 

(N 

CO 



CD 



i 

C 

o 
o 



> 
O 

OS 
O 



X 



Shaffique Adam and M. D. Stiles 
Center for Nanoscale Science and Technology, National Institute of 
Standards and Technology, Gaithersburg, Maryland 20899-6202, USA 
(Dated: September 8, 2010) 

Assuming diffusive carrier transport, and employing an effective medium theory, we calculate 
the temperature dependence of bilayer graphene conductivity due to Fermi surface broadening as a 
function of carrier density. We find that the temperature dependence of the conductivity depends 
strongly on the amount of disorder. In the regime relevant to most experiments, the conductivity 
is a function of T/T* , where T* is the characteristic temperature set by disorder. We demonstrate 
that experimental data taken from various groups collapse onto a theoretically predicted scaling 
function. 
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I. INTRODUCTION 

Monolayer and bilayer graphene are distinct electronic 
materials. Monolayer graphene is a sheet of carbon in a 
honeycomb lattice that is one atom thick, while bilayer 
graphene comprises two such sheets, with the first lat- 
tice 0.3 nm above the second. Since the first transport 
measurements^* 2 - in 2005, we have come a long way in 
understanding the basic transport mechanisms of carri- 
ers in these new carbon allotropes. (For recent reviews, 
see Refs.ii). 

A unique feature of both monolayer and bilayer 
graphene is that the density of carriers can be tuned con- 
tinuously by an external gate from electron-like carriers 
at positive doping to holes at negative doping. The be- 
havior at the crossover depends strongly on the amount 
of disorder. In the absence of any disorder and at zero 
temperature, there are no free carriers at precisely zero 
doping. However, ballistic transport through evanescent 
modes should give rise to a universal minimum quantum 
limited conductivity cr m m m both monolayer^ and bi- 
layer graphene^— The "ballistic regime" should hold so 
long as the disorder-limited mean-free path is larger than 
the distance between the contacts! 1 ^ 11 At finite temper- 
ature, the thermal smearing of the Fermi surface gives a 
density n(T) ~ T 2 for monolayer graphene. For ballistic 
transport in these monolayers, the conductivity a ~ yfnf 
for large n, so a(T) ~ T, 12 : 13 In the absence of disorder, 
cr(T) interpolates from the universal <r m j n to the linear 
in T regime following a function that depends only on 
T/Tp; (Tp is the Fermi temperature)^ 

Most experiments, however, are in the dirty or diffu- 
sive limit, which is characterized by a conductivity that 
is linear in density (i.e. a = ne/i c , with a mobility \x c 
that is independent of both temperature and carrier den- 
sit y i ), and the existence of a minimum conductivity 
plateau^ in cr(n), with cr m in = n lmB efj, c /V3. n lms is the 
root-mean-square fluctuation in carrier density induced 
by the disorder. In bilayer graphene, to our knowledge, 
all experiments are in the diffusive limit. 

The purpose of the current work is to calculate the 
temperature dependence of the minimum conductivity 



plateau in bilayer graphene. The temperature dependent 
conductivity of diffusive graphene monolayers is under- 
stood to depend largely on phonons^ but monolayer 
and bilayer graphene are distinct electronic materials and 
phonons are not expected to be important for bilayer 
graphene transport at the experimentally relevant tem- 
peratures^ 



II. THEORETICAL MODEL 

An important difference between monolayer and bi- 
layer graphene is the band structure near the Dirac point. 
Monolayer graphene has the conical band structure and a 
density of states that vanishes linearly at the Dirac point. 
Bilayer graphene has a constant density of states close to 
the Dirac point from a hyperbolic dispersion. The tight- 
binding description for bilayer graphen o 20 ' 21 results in a 
hyperbolic band dispersion 



E-p{n) = Vprn y/l + n/no — 1 



(1) 



that is completely specified by two parameters, vp w 
1.1 x 10 8 cm/s and n = v F m 2 /{h 2 Tr) « 2.3 x 10 12 cm" 2 
(where h — 2ttH is Planck's constant). For very small 
carrier density n <C fio> one can approximate bilayer 
graphene as having a parabolic dispersion, although most 
experiments typically approach carrier densities as large 
as 5 x 10 12 . The density of states for bilayer graphene is 
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where the parabolic approximation keeps only the first 
term. 

Understanding the temperature dependence of the con- 
ductivity minimum is complicated for two reasons. First, 
there is activation of both electron and hole carriers at 
finite temperature. Second, the disorder induces regions 
of inhomogeneous carrier density (i.e. puddles of elec- 
trons and holes). Moreover, tuning the carrier density 
with a gate changes the ratio between electron-puddles 
and hole-puddles, until at very high density there is only 
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FIG. 1: (Color online) Bilayer graphene mobility as a func- 
tion of back-gate voltage V g , normalized by the mobility at 
V g — 40 V. Solid lines use bilayer graphene's hyperbolic 
dispersion relation, while dashed lines are the parabolic ap- 
proximation valid only for low carrier density. Upper panel 
- long-ranged Coulomb impurities. From bottom to top: 
over-screened (parabolic), RPA (parabolic), Thomas- Fermi 
(parabolic), over-screened (hyperbolic), Thomas- Fermi (hy- 
perbolic). Lower panel: short-range (i.e. "delta-correlated" 
or "white noise") impurities. From bottom to top: RPA 
(parabolic), Thomas- Fermi (parabolic), Thomas- Fermi (hy- 
perbolic), unscreened (hyperbolic), unscreened (parabolic). 
See Ref. for definitions of the different approximations. 



a single type of carrier. The temperature dependence 
of the conductivity for bilayer graphene was studied in 
Ref. Ull using a coherent potential approximation. While 
this approach better captures the impurity scattering and 
electronic screening properties of graphene, it does not 
account for the puddle physics which is our main fo- 
cus. Reference [l6| modeled the temperature dependence 
of the Dirac point conductivity by assuming that the 
graphene samples comprised just two big "puddles" each 
with the same number of carriers. In the appropriate 
limits, our results agree with these previous works. Be- 
low we will provide a semi-analytic expression for the 
graphene conductivity by averaging over the random dis- 
tribution of puddles with different carrier densities. This 
result is valid throughout the crossover from the Dirac 
point (where fluctuations in carrier density dominate) 
to high density (where these fluctuations are irrelevant), 
both with and without the thermal activation of carriers. 

Given a microscopic model for the disorder, one can 
compute both \x c and n rms . Shown in Fig. [T] are re- 
sults for bilayer graphene mobility assuming both short- 
range and Coulomb disorder with different approxima- 
tions for the screening, and for both parabolic and hy- 
perbolic dispersion relations. As seen from the figure, 
generically, Coulomb impurities show a super-linear de- 
pendence on carrier density while short-range scattereres 



are sub-linear. Similar to monolayer graphene^^r— in- 
creasing the dielectric constant tends to decrease (in- 
crease) the scattering of electrons off long (short) range 
impurities, except in the over-screened and unscreened 
limits. All experiments to date find the mobility to be 
linear in gate voltage, so it is unclear what the dominant 
scattering mechanism in bilayer graphene is (see also dis- 
cussion in Ref. 1251 ) . Further experiments along the lines 
of Refs. iHl are needed. 

In what follows we take /i c and n rms to be parame- 
ters of the theory that can be determined directly from 
experiments: [i c can be obtained from low tempera- 
ture transport measurements and n rms from local probe 
measurements^ - — Lacking such microscopic measure- 
ments for the samples we compare with, we treat n rms 
as a fitting parameter, while taking [i c from experiment. 
As a consequence of this parameterization, the results re- 
ported here do not depend on the microscopic details of 
the impurity potential, provided this parameterization 
reasonably characterizes the properties of the impurity 
potential. Until more information about the important 
scattering centers is determined from experiment, all mi- 
croscopic models will require a similar number of param- 
eters such as the concentration of impurities ni mp and 
their typical distance d from the graphene sheet. Fur- 
ther, the results will disagree with experiment unless the 
choices give a constant mobility. 

A key assumption in this work is the applicability of 
Effective Medium Theory (EMT), which describes the 
bulk conductivity ctemt of an inhomogeneous medium 
by the integral equation^ 



dnP[nf (n) - aEMT = (>■ 



a(n) 



■ CTEMT 



(3) 



P[n] is the probability distribution of the carrier den- 
sity in the inhomogeneous medium - positive (negative) 
n corresponds to (electrons) holes, and cr(n) is the local 
conductivity of a small patch with a homogeneous car- 
rier density n. Ignoring the denominator, Eq. [3] gives 
cemt equal to the average conductivity. The denomi- 
nator weights the integral to cancel the build-up of any 
internal electric fields. The EMT description has been 
shown to work well whenever the transport is semiclas- 
sical and quantum corrections and any additional resis- 
tance caused by the p — n interfaces between the electron 
and hole puddles can be ignored £2r£2. It is assumed that 
the band structure is not altered by the disorder, which is 
to be expected for the experimentally relevant disorder 
concentrations*^ Since we are concerned with diffusive 
transport in the dirty limit, we expect that the EMT 
results hold for bilayer graphene. 



III. RESULTS 

To solve Eq.|3]we make the additional assumption that 
the distribution function P[n, n g \ is Gaussian centered at 
n g , (i.e. the field effect carrier density induced by the 
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FIG. 2: (Color online) Conductivity vs. gate voltage for clean 
and dirty graphene bilayers calculated from Eq. [3] Solid 
curves use the hyperbolic dispersion relation while dashed 
lines (only distinguishable at high temperature) show the 
parabolic approximation. Choice of parameters were based 
on experiments of Ref. [lU (clean) and Ref. H3 (dirty). Left 
panel: fi c = 6, 750 cm 2 /Vs, n rms — Ax 10 11 cm~ 2 and (from 
bottom to top) T = 20 K, 100 K, 180 K and 260 K. Right 
panel: fi c = 1,100 cm 2 /Vs, n rms = 1.25 x 10 12 cm -2 and 
(from bottom to top) T = 12 K, 105 K, 171 K and 290 K. 
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FIG. 3: (Color online) Minimum conductivity as a function of 
temperature for linear dispersion (upper curve) and parabolic 
dispersion (lower curve) graphene. Dashed lines show the 
high temperature asymptotes <7 m i n —¥ ire(i c T 2 / (3h 2 Vp) for 
linear and er m m — > me/x c 41n 2T/(nh 2 ) parabolic cases. Solid 
(red) line shows the hyperbolic result for n rms = 10 12 cm~ 2 . 
Also shown is that the hyperbolic result extrapolates from the 
parabolic theory at large a = rn 2 Vp/(h 2 Tm TlnB ) becoming sim- 
ilar to the linear dispersion for small a. (Red squares show 
results for a = 100 and red circles are for a — 0.01; here we 
ignore the contribution from higher bands). 



FIG. 4: (Color online) Same results as in Fig. [3] showing com- 
parison with experimental data from several groups. Inset 
shows the unsealed experimental data, while the main panel 
shows that the data collapses onto the theoretical curve with 
one scaling parameter (n rms ), where for each of these sam- 
ples, we also use the value of mobility reported by the authors 
and obtained from a separate low temperature measurement. 
Green triangles show suspended bilayer data from Ref. |401 us - 
ing /x c = 1.4 m 2 /Vs and T* = 36 K. Orange squares (Ref.[3J) 
and diamonds (Ref. lit! ) are bilayers on a Si02 substrate with 
fx c = 0.11 m 2 /Vs, T* = 530 K and fi c = 0.045 m 2 /Vs and 
T* = 290 K. Cyan circles show the four data points of Ref. ITsl. 
with /i c = 0.675 m 2 /Vs, T* = 80 K, which are off-scale in the 
main panel. 



back gate that is proportional to V g ), with width n rms . 
(This assumption is justified both theoretically^— and 
empirically 26 ). Our results are shown in Fig. [5J where 
as discussed earlier, the temperature dependence comes 
from the smearing of the Fermi surface. 

At first glance, it is not obvious that the results for 
clean bilayer graphene (left panel of Fig. [2]) and dirty 
bilayer graphene (right panel) are closely related. How- 
ever, if we consider scaling the conductivity as ctemt = 
CEMT/(^rmse)Ltc), scaling temperature as t = T/T*, 
where we define k^T* = Ep(n = n rms ), and scaling car- 
rier density as z = rt/n rms , we find that for both the 
linear band dispersion (n uq) and the parabolic band 
dispersion (n <C no), the scaled functions <5emt(z, t) each 
follow a universal curve. This is illustrated in Fig. [3] 
where we show the temperature dependence of the mini- 
mum conductivity. The results for the hyperbolic disper- 
sion (which is the correct approximation at experimen- 
tally relevant carrier densities) , depends on an additional 
parameter a — no/nrms)^ 

The scaling function for the hyperbolic dispersion ex- 
trapolates from the parabolic theory at large a becoming 
similar to the linear result for small a. For the experi- 
mentally relevant regime a ~ 1 the hyperbolic result de- 
pends only weakly on a and is indistinguishable from the 
parabolic result for T < 0.5 T*. 
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of a universal scaling form, we consider the conductivity 
for a parabolic band dispersion. Using the scaled vari- 
ables defined above, we can manipulate Eq. (j3|) into the 
dimensionless form 



dzexp [— z 2 /2] cosh[z z 

o 9 #M 



= 0, (4) 



n r 



FIG. 5: (Color online) Bilayer layer graphene conductivity 
as a function of temperature and carrier density for T/T* — 
0,0.5, 1, 1.5, and 2. Inset shows a close-up of the zero tem- 
perature minimum conductivity (which is the same for both 
monolayer and bilayer graphene). The dashed horizontal line 
shows the result for <7 m i n , while the other dashed line is the 
high-density transport regime. The solid line (Eq. U) captures 
the full crossover from the regime where the conductivity is 
dominated by the disorder induced carrier density fluctua- 
tions, to the semiclassical Boltzmann transport regime. 



This analysis suggests that <7 m i n (T)/ (e/x c ), which can 
be taken directly from experiment, is not a function of 
/j, c , but only n rms . We take results from a set of exper- 
iments in very different regimes (see the inset of Fig. [4j 



where z g = n g /n vms and we have written the local con- 
ductivity as a(n,T) — n Tms efi c H(z,t). Below we cal- 
culate the dimensionless function H{z,t) assuming ther- 
mally activated carrier transport with constant n rms and 
fi c and explicitly show that it depends only on scaled 
variables z = n/n Tms and t = T/T* . With the analytical 
results for H(z,t) discussed below, this implicit equation 
can be solved either perturbatively or by numerical inte- 
gration to give (Temt- The results of this calculation are 
shown in Fig. [5] 

To proceed, we calculate the function H(z,t). For ther- 
mal activation of carriers, the chemical potential \i is de- 
termined by solving for n g = n e — n^r^ where 

/■oo 

n e (T) = / dE D(E)f(E,v,k B T), 
Jo 

n h (T) = f dE D(E)[l-f(E,fx,k B T)], (5) 

J — oo 

where f(E,fj,,k B T) is the Fermi-Dirac function and k B 
is the Boltzmann constant. For T = 0, only major- 
ity carriers are present, while for T — > oo, activated 
carriers of both types are present in equal number. 
Within the parabolic approximation, we find n e (/o = 
,(T/T F ) In [1 + exp(Tfi/k B T)} and fi = E F . Using 



and choose n rms to fix the value of CT min (T)/(n rms e/Xc) at T ) = (" e + n hWc we obtain 

H(z,t) = z + 2iln [l 



T = 0. Then using k B T* (n rms ) = ^(rirms) to scale the 
temperature, all of the results lie on top of the theoreti- 
cal curve computed using the hyperbolic dispersion, see 
Fig. [4] The theoretical curve with which they agree is dis- 
tinct from similar curves calculated for a linear dispersion 
and for the purely parabolic dispersion at high T/T*. We 
note that the scaling function is more complicated than 
a line. The calculation reproduces not only the initial 
slope as a function of temperature, but the crossover to 
higher temperature behavior. For the parabolic disper- 
sion, which agrees at low temperatures, the conductivity 
extrapolates from <7 m i n (T — > 0)/(n rms e/x c ) s» at 
low temperature to er m i n (t 3> l)/(n rms e// c ) ~ (21n2)i at 
high temperature, with a crossover temperature scale of 
T m T*/2. In the future, it should be possible to fur- 
ther test this agreement by measuring n Yma experimen- 
tally^^ 

One feature of Fig. [2] and Fig. [4] is that for most of the 
experimentally relevant regime, the temperature depen- 
dence of the conductivity calculated using the parabolic 
approximation provides an adequate solution. This limit 
has been treated in contemporaneous wor k 41 i 42 treating 
this problem with different approximations and reaching 
similar conclusions. To better understand the emergence 



,-z/t 



(0) 



This demonstrates that Eq. 2] depends only on the scaled 
variables, guaranteeing that ctemt is a function only of 
T/T* and n g /n vms as shown in Fig. [SJ 

A similar analysis can be done for the hyperbolic dis- 
persion. We find 



H(z, t, a) = 



C + 2 



(7) 



where g(z,a) = T*/Tp, £(z, a) = — 1 + y/T+z/a, and 
the scaled chemical potential y — fi/E-p is given by 



V= 2 



2 + e - 2atg) 2 (U 2 (-e-y^) - Li 2 (-e+^))l ,(8) 



where Li2(z) = f dt t^ 1 ln(l — t) is the dilogarithm func- 
tion. Only for a 3> 1 and a <C 1 does H(z,t,a) become 
independent of a = uq /n rms giving the universal scaling 
forms for linear and parabolic dispersions, respectively. 
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IV. CONCLUSION 

In summary, we have developed an effective medium 
theory that captures the gate voltage and temperature 
dependence of the conductivity for bilayer graphene. The 
theory depends on two parameters: n rms that sets the 
scale of the disorder, and /i c the carrier mobility. These 
could be computed a priori by assuming a microscopic 
model for the disorder potential and its coupling to the 
carriers in graphene. Alternatively, one could use an 
empirical approach where one uses experimental data at 
T = to determine the parameters and use the theory 
to predict the temperature dependence. 

Our main finding is that experimental data taken from 



various groups collapse onto our calculated scaling func- 
tion where the disorder sets the scale of the temperature 
dependence of the conductivity. This further suggests 
that even some suspended bilayer samples are still the 
the diffusive (rather than ballistic) transport regime. 
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